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Abstract 

Ambiguities in the functional-integral solution of the stochastic differential equation (SDE) 
arising due to the definition on the functional Jacobi determinant and the white-in-time limit 
in the noise are analyzed and two forms of the de Dominicis-Janssen dynamic action proposed 
corresponding to the Ito and Stratonovich interpretations of the SDE. 

Introduction 



Since the invention of the functional-integral form [T] of the Martin-Siggia-Rose approach to solution of 
stochastic problems [2] ambiguities in this formulation have been discussed in the literature O [H O [6] . 
. They have two sources: first, a change of variables brings about a Jacobi determinant which is not 
C ! well defined. Second, the mathematical inconsistency of the Langevin equation with the white noise 
' leads to different definitions of the SDE that show in the functional-integral solution. The aim of this 
^ ■ report is to analyze both cases in unified manner. Two different ways to resolve the ambiguities in the 
^ \ functional integral are proposed in the terms of dynamic actions in the Ito and Stratonovich forms. 

! Dynamic field theory with additive random field 

QQ ■ Consider the Langevin equation 

^ ■ r) 

-^=V{^) + f:=-K^ + U{^)+f, (1) 



where / is a Gaussian random field with zero mean and the white-in-time correlation function 

(/(t, x)/(t', x')) = D{x - x') = 5{t - t')D{^ - x') , (2) 



^ ! where the shorthand notation x = (t, x) has been used. 



In ([T]) is usually a second-order differential operator in space and U{x, ip) a nonlinear local in 
time functional of ip. The paradigmatic example is model A of critical dynamics defined by 



^ = _r i^-vv + + g ^'J + / , (3) 

and the correlation function D{x — x') = 2T6{x — x'). 

The Langevin equation with white-in-time noise / is mathematically inconsistent, because the 
time integral of the noise J fdt is a Wiener process which not differentiable anywhere. To arrive at 
this limit, use the set of correlation functions constituting a 6 sequence in time, i.e. 

{fit, x)f{t\ x')) = D{t, x; t\ x') — ^ 8{t - t')D{x, x') . (4) 
This treatment gives rise to the solution of the SDE ([T]) in the Stratonovich sense [S] . 
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To cast the problem ([T]) and (jlj) into a functional-integral form, use standard tricks: 



Vif {5 [-dtif + V{ip) + /] I det M|) e 



= yp(^ Jvip (^\detM\ e^[-9*^+^(^)+/l ^ e^-^ . (5) 
Representation the Jacobi determinant in ([5|) in the loop expansion of eT--!-^^ yields 



det M = det 



5ip{x') 



det {dt + K) e 



-A{0)U' 



(6) 



in which the diagonal value A(0) of the response function of is a parameter of the model. The 
response function is the retarded Green function of the free-field equation, whose value at t = t' is 
not determined. The two popular choices are A(0) = [S] and A(0) = ^(5(x — x') [4J. Contrary to 
claims in the literature [3l [5], the determinant ([6|) has nothing to do with the stochastic nature of 
the Langevin equation ([1]): it appears in the functional-integral representation of the solution of the 
corresponding deterministic equation as well. 

Thus, we arrive at the De Dominicis-Janssen action 



Dominicis-Janssen action 
= -A(0)[/' + In P[f] + ^ [-dt^ -Kip + U{ip) + f] , 
density function of the field f . The effect of the de 



(7) 



where P[f] is the probability density function of the field /. The effect of the determinant ([6]) is to 
cancel superfluous graphs brought about by the usual Feynman rules for the dynamic action ([7|). For 
instance, the iterative solution of ([3]) for the field may be expressed as the tree-graph expansion 



+ 




+ 



(8) 



Here, the directed line corresponds to the propagator A, the cross to the field / and the full dot to the 
vertex factor |9j of the dynamic action d?]). However, according to the Feynman rules for ([7|), in the 
linear order in A, this solution should contain also a graph with a closed propagator loop and graph 
brought about by the determinant factor (here, the source field / is fixed and lnP[/] = 0), i.e. 




The black square stands for the vertex factor brought about by the determinant term — A(0)C/' = 
A(0)^Ar(/3^. The coefficient A(0) is the value of the closed propagator loop of the third graph in the 
right side of therefore the third and fourth graphs in ^ cancel each other thus restoring the 
iterative solution ([8]) at this order. 

Although the functional integral gives rise to the perturbative solution independent of the value 
of A(0), the form of the dynamic action is important for its nonperturbative calculation. The choice 
of det M affects the explicit form of the integrand of the functional integral in a nontrivial way. The 
subsequent implicit dependence of the functional integral on the parameter A(0) should be dealt with 
in nonperturbative approaches such as numerical evaluation or instanton calculus |10l [TT] to arrive at 
generating function independent of A(0). To be on the safe side, for purposes of any nonperturbative 
calculation, the Jacobian should be taken in the field-independent form det M = det (dt + K) , which 
corresponds to the choice A(0) = 0. 
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Dynamic field theory with multiplicative random field 

A different picture emerges for a SDE with a multiplicative noise. Consider the simplest case 

^ = -K^ + f^, (10) 

where X is a time-independent operator acting on the field and / a Gaussian random field with 
zero mean. Difficulties in the interpretation of the equation (jlOp arise, when the correlation function 
of the field / is local in time as in ([2]), which is limiting case of the sequence (jl]). 
The iterative solution of the SDE ([10]) may expressed as the series 



= Ax + A/Ax + A/A/Ax +... (11) 

where x is the initial condition of the solution Ax = / dx.'A{t, x— x')x(x') of the homogeneous equation 
(llOp where A is the (retarded) Green function of the same equation. Graphically, the solution (llip is 
a sum of chains of oriented lines corresponding to retarded propagators 



■ \ ■ o + . \ . \ 



^=^-0 + . f -.O + ' «i ' O + ... (12) 



where the circle stands for the initial condition x of the homogeneous equation (jlOp . the wavy line 
corresponds to the random field / and the full dot represents the vertex factor brought about by the 
last term of (jlOp . The functional-integral representation ([5]) for the solution of this problem gives rise 
to the dynamic action 

S[^, if, f] = -A(0)/ + In P[f] + f i-dtif -Kf + fif) (13) 

again with a term brought about by the determinant factor det M = det {dt + K) exp [— A(0)/] which 
has the same effect of cancelation of graphs as in the case of additive external field. The perturbation- 
theory expression for the field f at the order depicted in (|12p is 



^ = — O + ' « « ' O + 



. O -I- . \ 




2 ^2 

which again contains closed loops of the propagator. Here, the black square stands for the — A(0) 
vertex factor. Thus, the sum of the closed propagator loop graphs and the determinant-generated 
graphs in the right side of (114p vanishes. Note that the effect of the Jacobi determinant to remove 
"superfiuous" graphs with closed propagator loops is independent of the nature of the source field /. 

The average of the series (llip over the probability density of / gives the solution of the SDE (llOp 
for the field (p. For the Gaussian distribution with zero mean and the correlation function ([2]) this 
gives rise to the dynamic action 

S[ip, ^, /] = - A(0)/ - \fD-'f + f i-dtip -Kf + fip) . 

The analysis of the corresponding functional integral is more transparent in terms of the three fields 
if, ip and /. The perturbative solution of the SDE pUj) is given by Wick's theorem for the Gaussian 
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distribution of /, which graphicahy amounts to replacing any pair of / by the correlation function D 
depicted by an unoriented line in all possible ways. For instance, 

(^)=_0 + . <^ . O + ... (15) 



In the limit of white-in-time correlations an enormous truncation of the averaged iterative series (jlip 
takes place due to temporal 5 functions contracting the ends of chains of the retarded propagators. 
Only those terms, in which the correlation function is multiplied by a single retarded propagator do 
not vanish automatically. 

For instance, the one-loop graph in ()15p in case of white-in-time noise gives rise to an ambiguity, 
which is directly related to that in the interpretation of the SDK (jlOp . A straightforward substitution 
of the white-noise correlation function in this graph gives rise to the expression 



j dti j dxi J dx2 J (ix3 A{t — ti, x — xi) 



X A(0, Xi - X2)I?(xi, X2)A(ti, X2 - X3)X(X3) , 

where the value of the propagator at coinciding time arguments A(0,xi — X2) = 6{0)5{xi — X2) is 
again ambiguous, although here the spatial 6 function does not create a problem. 

With the use of the (5-sequence of correlation functions @ this ambiguity is readily resolved and 
gives rise to the expression 



^ y"c?ii y dxi J dx3A{t - ti,x - xi)Z?(xi,xi)A(ti,xi - X3)x(x3) , (16) 



with the coefficient i in front of the spatial 5 function. Formally, this result may be obtained by 
amending the definition of the propagator according to the rule A(0,x — x')) = ^(5(x — x'). Then, 
of course, it would seem appropriate to use the same rule in the Jacobi determinant factor. Within 
this choice the graphical expression in (jl6p appears excessive, because it hints to twice the number of 
integrations than actually is carried out. 

It is more convenient to adopt the convention A(0) = here as well. In this case the deterministic 
part of the dynamic action (|13p has to be amended to produce the correct nonvanishing result for 
graphs like (jl6p . This may be effected by adding to the dynamic action the term 

AS" = ^ (pD{0)ip ■.= ^ jdt jdx (p{t, x)D{x, x)ip{t, x) . (17) 

Graphically, this amounts to replacing the one-loop graph with the noise correlation function by a 
new vertex factor with the coefficient ^-D(O): 



It may be readily seen that the amendment (I17p to the dynamic action corresponds to the reorganiza- 
tion of the SDE (jlOp . when the Stratonovich equation is replaced by the Ito equation with the same 
solution [8j. Thus, we arrive at the dynamic action in the Ito form 



S'[^,^,f] = -^fD'^f + ^ 



-dtip-Kip+^DiO)^ + fip 



A(0) = 0. (18) 
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On the other hand, if we adopt the choice A(0, x — x') = ^(5(x — x'), then it is natural to make the same 
choice in the determinant factor and include it in the dynamic action. In this case, the determinant 
term in the dynamic action leads to cancelation of terms containing diagonal values of the propa- 
gator A(0) J / dx J dx.'D{-K — x')A(O) = jS'^{0)TV f dxZ)(x) but leaves intact terms containing a 
similar expression but with diagonal in time values of propagators contracted with different spatial 
arguments: / / dx / (ix'A(0,x - x')-D(x - x')A(0,x' - x) = ld{0)TVD{0) . To exclude the latter 
from perturbation theory, the term AS' = — | 6{0)TVD{0) must be added to the determinant part of 
the dynamic action. In this case we arrive at the dynamic action in the Stratonovich form 



5"b,(^,/] = -i<5(0) 



^^(0) + / 



\fD-'f + ^ [-dtif -K^ + fip] . (19) 



In case of multiplicative noise the choice of the diagonal value of the propagator affects the dynamic 
action in a fashion analogous to the relation between Ito and Stratonovich SDE's. It should be 
emphasized that both dynamic actions (jlSp and (jl9p give rise to the same correlation and response 
functions which correspond to the Stratonovich solution of the SDE (llOh . With the use of similar 
arguments, two dynamic actions may be constructed for the solution of the SDE in the Ito from as 
well. Different dynamic actions for the functional-integral solution of a SDE represent ambiguities in 
the construction of the functional integral, but not in the SDE itself. 

The analogy between the compilation of dynamic actions and SDE holds also in the case, when 
the random source field is multiplied by a polynomial function(al) of the field (p. Consider the SDE 

dt ^^2-'^ 
as an example. The dynamic action with the coloured noise is 

1 -1 . . ~ / . .1.2 



S\^, (f>, f] = -A(0)/<^ --fD f + ^l -dtv -Kip + -f^'] , (20) 



where the first term in the right side is brought about by the determinant factor. Instead of (jlSp the 
first terms of the solution are (take the d sequence of the correlation functions of /) 



{<^)=— o + ■ ^ ^^^^ 

The analytic expression of the second term in the white-noise limit assumes the form 



O ~ 2 J J ^'^^ J '^'^^ J ^'^^ J ^'^^ ~ ^i)X — xi) 



X L>(xi,Xi)A(fi,Xi - X2)x(x2)A(ti,Xi - X3)x(x3)A(fi, Xi - X4)x(x4) . 

The same result is brought about by the term AS = ^(pD{0)ip^ added to the dynamic action f2U\) . 
The form of this term corresponds to the transition from the Stratonovich SDE to the Ito SDE. 
Graphically, the one- loop white-noise subgraph shrinks to a dot: 




In case of Stratonovich convention, terms with propagator loops contracted with different spatial 
arguments are accounted for by the term AS" = — |(5(0)D(0) J dt f d'Kip'^{x) , which should then be 
added to the dynamic action (PU|) . 
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Iterative solution of the stochastic differential equation 

Consider the SDE with a multiphcative white-in-time noise in the mathematically consistent form W 
of the form 

dip = -Kipdt + b{ip)dW . (22) 

Here, dW is the increment of a Wiener process for each value of the spatial coordinate x, the deter- 
ministic part is linear in ip for simplicity whereas the coefficient function 6(x, cp) of the stochastic part 
is a polynomial time-local functional of (p. 

Let us construct a solution for the SDE (f22]) by iterations. To this end, introduce first spatial 
Fourier transform to obtain 

t t 

ifit, k) = v^(0, k) - y" Vit', k)dt' + b{ip, q)dW{t', k - q) , (23) 



where K is now a constant and b{ip, q) is the Fourier transform of the coefficient function b{ip) of (|22p 
and dW{t, k) is the that of the increment of the Wiener process. In ()23p the function b{(p) is expressed 
as a functional of the field ip in the Fourier representation. Construct then the iterative solution for 
the Laplace transform 0(s, k) of the field ip{t, k). The first two steps of the iteration process yield 





Here, (p{0, k) is the Fourier transform of the initial value of the problem. The iterative solution of (I23p 
is readily expressed in terms of a tree-graph expansion of the type (|12p , when the coefficient function 
b{(p) is a polynomial function of tp. 

The ambiguity in the solution of (j22p shows, when averages of functions of Lp are calculated. The 
combinatorial rule for the Wiener process is given by the usual Wick theorem. The result of calculation 
of time integrals of pairwise products of the Wiener process in case of Ito stochastic integral is given 
by the correlation formula ^ 



G{u)dW{u) J H{v)dW{v)j = J{G{u)H{u))du, 

to to to 

where G{t) and H{t) are continuous nonanticipating functions and the brackets denote average over 
the probability distribution of the Wiener process. All functions of time emerging in the iterative 
solution of the SDE (j22p are nonanticipating. 

The practical rule of calculation may be stated as follows [8]: rewrite the increment of the Wiener 
process as dW(t) S,(t)dt , where ^{t) is (5-correlated in time, i.e. (C(t)'^(t')) = 6{t — t') , and calculate 
the averages in terms of ^. With the use of this rule, integrals appear in which the 6 function of 
correlations of ^ is integrated in such a way that one of the limits of integration coincides with one of 
the arguments of the 6 function. With the use of Ito stochastic integral, ambiguities arising in these 
case are resolved as 

t2 t2 

fmt - h) = fih) , / f{t)S{t - t2) = . (24) 



tl tl 
In case of the Stratonovich stochastic integral these rules are replaced by 



t2 t2 

fmt - h) = Ifih) , f fmt - = lf{t2) . (25) 



tl 
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In the tree-graph iterative solution these situations arise solely due to the presence of the temporal 
step function in the propagator. The Ito rules ()24p exactly correspond to the choice 6{0) = 0, whereas 
the Stratonovich rules ([251) correspond to the choice 0(0) = | for the value of the step function at the 
origin. 

Thus, we see that averages of the iterative solution of the SDE (j22p in the Stratonovich sense give 
rise to the same expression as the perturbation expansion on the basis of the dynamic action in the 
Stratonovich form: 

S"[^, ^, f] = -^7d?-^ b'f - \fD-^f + ^ [-dt^ -Kip + fb{ip)] , A(0, X - x') = ^<5(x - x') . 
Here, the determinant contribution in an expanded form is 

'-7m + '-b'f='- [dt /dx/dx'|^Z?(x-x')?^ + - [dt [dJ-pfl4fit,x). 
8 2 ^ 8j J J 6ip{t,x') ^ ' 6^{t,x) 2 J J Sip{t,x)^^' ' 

It has been shown above that the same solution may be obtained in the perturbation expansion based 
on the dynamic action in the Ito form. A generalization of that argument to the case of arbitrary 
monomial shows that for the case of polynomial b{^p) the additional term needed to obtain the Ito 
form of the dynamic action corresponding to the SDE (I22p is 

AS=\^b'{^)Db{^) = \ jdt jdxjdx!^{t,x)^^^^D{x-x!)b{x!,^). (26) 

Here, Z)(x — x') is the correlation function of the Wiener processes as functions of the coordinate, i.e. 
(VF(t,x)T^(t,x')) = tD{-x. — x') instead of the relation {W{t)W{t)) = f for a single Wiener process. 

Thus, the iterative solution of the SDE (|22p in the Stratonovich sense may be constructed with 
use of the Ito dynamic action of the form 



-dtv -K^ + \ b'{ip)Db{^) + fb{^) 



A(0,x -x') =0, 



where we have used the condensed notation (1261). 



Dynamic action in the Ito form and the Stratonovich form 



The analysis of preceding sections allows to write down the following prescriptions for the construction 
of the solution of the following Stratonovich SDE 



dt^ = -K^ + U{ip) + fb{ip) 



(27) 



in which the deterministic part is written in a generic form with the nonlinear term U{(p). 

The solution of (|27p is brought about by the functional-integral representation either with the 
dynamic action in the Ito form with A(0,x — x') = 



S'[^,^J] = -^fD-'f + if 



-dtif -Kip + U{^) + ^ b'{^)Db{^) + fbiif) 



or with the dynamic action in the Stratonovich form with A(0,x — x') = i(5(x — x') 



u'{^) + \7Dy+b\^)f 



-JD-^f + if [-dtf -K^ + U{ip) + fb{ip)] . 



From these representations the Ito is both more convenient in practical calculations in perturbation 
theory and more reliable for use in non-perturbative evaluation of the functional integral. 

Analogous representations for the action functional in stochastic dynamics defined by the Fokker- 
Planck equation have been obtained also within the operator approach of quantum field theory [12] 
and interpolation construction of the path integral 
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Conclusion 

Ambiguities appearing in the functional-integral solution of stochastic differential equations have been 
analyzed. The representation of the generating function of the differential equation as a functional 
integral gives rise to an ambiguity, which is completely independent of the statistics of the source 
field: this ambiguity in the choice of the initial value of the propagator appears even in the case of 
deterministic equation. The other ambiguity is due to the mathematical inconsistency of the SDE 
with white-in-time noise, whose resolution gives rise to interpretations in the Ito sense and in the 
Stratonovich sense. It has been shown that the functional-integral solution of both SDEs may be 
constructed with the aid of two dynamic actions, which differ by the choice of the initial value of the 
propagator and the deterministic part of the dynamic action. The difference of the deterministic parts 
of the two dynamic actions corresponds to the connection between Ito and Stratonovich SDEs having 
the same solution. 
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